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Abstract 

We report on results of computer simulations of spherical resonant-mass 
gravitational wave antennas interacting with high-frequency radiation from 
astronomical sources. The antennas were simulated with three-mode induc- 
tive transducers placed on the faces of a truncated icosahedron. Overall, 
the spheres were modeled with a sensitivity of about three times the standard 
quantum limit. The gravitational radiation data used was generated by three- 
dimensional numerical computer models of inspiraling and coalescing binary 
neutron stars and of the dynamical bar-mode instability of a rapidly rotating 
star. We calculated energy signal-to-noise ratios for aluminum spheres of 
different sizes cooled to 50 mK. We find that by using technology that could 
be available in the next several years, spherical antennas can detect coalescing 
binaries out to slightly over 15 Mpc, the lower limit on the distance required 
for one event per year. For the rapidly rotating star, we find, for a particular 
choice of the radius at centrifugal hangup, spheres are sensitive out to about 
2 Mpc. The event rate is estimated to be about 1 every 10 years at this 
distance. 

PACS number (s): 04.80.Nn, 95.55.Ym, 04.30.Db 
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I. INTRODUCTION 



The experimental effort to directly detect the space-time ripples known as gravitational 
waves has been going on for 35 years, beginning with Weber's pioneering work in the early 
1960s Since then, two main experimental approaches have evolved: cryogenic resonant- 
mass detectors [Q-^ and laser interferometers The sensitivity of both techniques is 
sufficient that unambiguous detection of gravitational waves is expected soon, perhaps within 
the next ten years. The work of Taylor and Hulse [^], showing the orbital decay of binary 
pulsar PSR 1913-16 is in agreement with general relativity's prediction for gravitational 
wave emission, has added to the anticipation of the first direct, confirmed detection. It 
appears possible that a new generation of advanced resonant-mass detectors could operate 
concurrently with interferometers already under construction. At this time, understanding 
possible sources of gravitational waves and which experimental technique is best suited 
to which source of radiation takes on greater importance. We endeavor to clarify this 
by numerically computing energy signal-to-noise ratios for resonant-mass detectors and 
interferometers interacting with two possible sources of gravitational waves. 

The best understood source of detectable gravitational waves is from inspiraling and co- 
alescing binary neutron stars ^ . The Laser Interferometer Gravitational-wave Observatory 
(LIGO) has been designed and optimized to detect these events at a distance of 200 Mpc 
after significant interferometer improvements [^]. To accomplish this, it has been devised 
to be most sensitive at as low a frequency as possible (~ 200 Hz) where the waveform from 
binary neutron stars is stronger. However, the waveform at 200 Hz is due almost solely to 
the inspiral phase of the binary neutron star evolution and contains virtually no information 
about the coalescence. As the inspiral is determined by point-mass dynamics, the equation 
of state for nuclear matter (i.e., neutron stars) will affect only the coalescence waveform. Co- 
alescence also occurs when the gravitational field between the neutron stars is strongest, so 
the effects of general relativity will be more important than during the inspiral. To measure 
these effects, it will be necessary to monitor the higher-frequency waves from coalescence in 
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addition to those at lower frequency from inspiral. 

Resonant-mass gravitational wave detectors have been in use for longer than the inter- 
ferometric detectors. Resonant-mass antennas with bar geometries have been taking data 
and been continually improved since their inception. The use of spherical geometry as 
an improvement over bars was first suggested by Forward in 1971 |T^, and Wagoner and 



Paik later showed that at equal frequencies spheres have an advantage over bars in energy 
cross section ||ll|. Recent re-examination of spherical detectors ||12|- |T5|] has generated in- 



ternational interest in constructing one or more massive spherical antennas incorporating 
advances in transducer technology. The possibility of building such an antenna to operate 



concurrently with the first LIGO interferometers appears good |T^. The sensitive frequen- 
cies for a sphere are higher than those for the first LIGO interferometers, spanning about 
750 Hz to 2700 Hz in the lowest mode, and therefore they are well suited to complement 
interferometer experiments at high frequencies. One such advanced resonant-mass detector 
can be more sensitive than the first LIGO interferometers within a bandwidth of around 
100 Hz centered at the quadrupole resonance of the sphere and would have a sensitivity 
within that band comparable to what will be achieved by the advanced LIGO interferom- 
eters. The sphere's sensitivity is independent of source direction, unlike LIGO. Spherical 
antennas can also provide direction and polarization information more easily than LIGO [l^ 



and, properly instrumented, could detect any scalar gravitational radiation that might be 
present |lTO| , pJ^ . In this paper, when we refer to "high-frequency" gravitational radiation, we 
mean those signals that include significant strength above 750 Hz. This frequency is where 
the first LIGO interferometers' sensitivity begins to weaken from photon shot noise in the 
lasers and the spherical resonant-mass detector's sensitivity becomes important. 

We have looked at the question of whether a spherical detector, or in particular a Trun- 
cated Icosahedral Gravitational- wave Antenna (TIGA) as described by Merkowitz and John- 
son |T^, is capable of observing high-frequency events. Specifically, the coalescence of a 
binary neutron star system and the dynamical bar-mode instability of a single, rapidly 
rotating star were examined as possible astronomical phenomenon that could produce high- 



frequency gravitational radiation. Waveforms for these events, generated with computer 
simulations by Centrella's group at Drexel University, were used as input into a mathe- 



matical model of a 50 mK spherical detector with three-mode inductive transducers []T7 
The energy signal-to-noise ratios obtained from this model help determine how TIGAs and 
interferometer experiments can best complement one another. 



Coccia and Fafone |I8| have also looked at energy signal-to-noise ratios from astronom- 
ical events in spherical detectors. Our work and theirs are complementary. They looked 
solely at inspiraling binary neutron stars as sources, leaving out the coalescence phase as 
well as any other high-frequency events. Since the inspiral can be modeled accurately by 
point-mass dynamics, they used an analytical expression for the waveform. We found it nec- 
essary to use numerical data from computer models to simulate the coalescence. By limiting 
themselves to the inspiral phase, Coccia and Fafone were unable to accurately predict energy 
signal-to-noise ratios for higher mass neutron stars or black holes. For some sphere sizes 
and compositions their simulation does not produce results for 1.4 Mq neutron stars, the 
observed mass for all known neutron stars in binaries . However, their method was able to 
show that spherical antennas can determine the chirp mass 

]TP|), the orbital inclination and the distance to the source, a result our method did not 
produce. Both techniques provide useful information that are unobtainable by the other. 

In Section |I|, we describe the method used for the energy signal-to-noise ratio calcula- 
tions, how the code that produced the results was written, and what parameters for the 
spherical antenna we used. In Section |IT|, we discuss the signal waveforms we used as inputs 
into the model developed in Section H and present the results of the calculations. Finally, 
in Section |^ we present our conclusions and discuss ideas for further work. 

II. METHOD 

To calculate the energy signal-to-noise ratio (SNR) per unit bandwidth of the TIGA, 
we followed the method of Price who showed how to calculate the SNR for a bar 
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antenna that uses an optimal filter to process the data. Stevenson |^ has shown that 
for six identical transducers in the TIGA geometry and identical Qs for all five quadrupole 
modes of the sphere, the SNR of a spherical antenna is identical to that of a bar antenna 
instrumented with one of those six transducers. The equivalent bar has an effective mass as 
seen by the transducer of 
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mes = -x{^T^R^)p, (1) 

where R and p are the sphere radius and density. The dimensionless parameter x comes 
from the radial driving point admittance matrix of the sphere at the quadrupole frequency. 
For an aluminum sphere with a Poisson ratio of 0.33, x = 0.301 [T^. The factor of 5/6 in 
Eq. (|l]) accounts for the multiple sphere modes and transducers . The SNR for the TIGA 



is the same as that for the equivalent bar provided one equates the energies deposited by 
the gravitational wave in the two antennas. We can now calculate the SNR for the simpler 
case of a bar, while retaining all the information available from a sphere. 
The energy deposited in the sphere is calculated from 

E = TJl, (2) 

where S is the energy cross section of the sphere [0, 

(3) 

Jo 

Here Vg is the extensional sound speed of the sphere material, /o is the quadrupole frequency, 
and n is a dimensionless constant that accounts for antenna geometry and mode. It has the 
value 0.215 |13| for a sphere in the lowest quadrupole mode and 0.585 in the first excited 



quadrupole mode ||13|JT5[| . Throughout G and c are Newton's gravitational constant and the 



speed of light. The total energy flux is |22 



T^ = '^^^^h{u)\\ (4) 

G lOTT 

where lo is the angular frequency of the gravitational radiation and \h[LS) \ is the magnitude of 
the frequency-domain amplitude of the gravitational wave. Thus the total energy deposited 
in the equivalent bar is 



E = lP^U\hM\'. (5) 
^ Jo 



We define an effective force that acts on the equivalent bar [|13 



/eff(t^) = -mcs^^'^h{uj)ies, (6) 

where the relationship between fgfr and R is determined as follows. For an impulsive force, 
the energy is deposited as kinetic energy in the antenna. The energy E after the impulse is 
given by 

Combining Eqs. (^Uh gives 



'n^. (8) 



Then using Eq. (|l|) and the relationship between ujq,Vs, and R p3[ for each quadrupole mode 



an ies of 0.337i? in the lowest quadrupole mode and 0.109i? in the first excited quadrupole 
mode are calculated. Using Eq. (|^) as the definition of a force on the equivalent bar, the 
method of Price can be followed exactly. 

The transducer we assumed was a three-mode inductive transducer. A three-mode trans- 
ducer is necessary, rather than the standard two-mode system, to get higher bandwidths, 
which are required to reach sensitivities near the standard quantum limit. Higher bandwidth 
reduces the requirement on the Q of the sphere and transducer. Higher bandwidth is also 
useful to cover a larger spectrum of frequencies and reduce the need for additional antennas. 

Assuming that a template of the gravitational waveform is available, optimal filtering 
can be used on the output signal of the transducer. The optimal filter produces the highest 
SNR possible and has the form 



where u{uj) is the velocity signal of the antenna effective mass and Sn{uj) is the total velocity 
noise spectral density, both referred to the input of the optimal filter. The parameter to is 
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the time at which the SNR is to be optimized. To calculate u{lj) and Sn{uj), it is necessary to 
solve the equations of motion for the antenna coupled to the three-mode resonant transducer. 
They have the form: 

JUrriefflil = tl (Mint - Ml) H Ml, 

UJ uo 

jkint / N J^trans/ \ r /in\ 

JCUmintMint = (Mint " Mi) (^2 " Mint) " J2, (10) 

J^trans / \ , 

JCc;mtransM2 = (M2 - Mint) + /2- 

UJ 

Here, meff is the effective mass of the antenna, mint is the mass of the intermediate resonator, 
and mtrans is the transducer mass; /ces is the effective spring constant of the antenna, feint is the 
spring constant that connects the antenna to the intermediate mass, and fetrans is the spring 
constant between the intermediate mass and the transducer mass. The spring constants are 
complex valued and include dissipation. The variables Mi,Mint, and Ui are, respectively, the 
velocities of the antenna surface at the transducer, of the intermediate mass, and of the 
transducer mass. Applied forces acting on the antenna surface and between mtrans and mint 
are denoted by /i and /2, respectively. 

Eliminating Mint from Eq. ([1^) allows the equations of motion to be written as 

Ui = yijfj, (11) 

with i and j taking the values 1 to 2. The energy SNR per unit bandwidth, cr(u;), then 
becomes [^| 



a{uj) = K{lu)u2 (12) 



From Eq. (pT]), cr(co') is found to be 



(13) 
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2 + 2fcBTRe(i/22) + 2Re(2/22^/J 



^(^) = C , C L. 12 , OL N , OP„^.. C ^ ' (14) 



assuming no force on the transducer, i.e., /2 = 0. Here the force fi is fes from Eq. (^) and 
T is the physical temperature of the sphere. The matrix yij{uj) is the admittance matrix of 



the antenna with transducer defined in Eq. (pi]). The four terms in the denominator are the 
individual parts of Sn{uj), the velocity noise. They are, respectively, the additive velocity 
noise, the force noise, the thermal noise, and the correlation noise. The spectral densities 
are defined as 



Sfiuj)^ / e-^--(/(t)/(t-r))rfr, 

J— oo 

/ + 00 
e~-'^"(n(t)M(t-r))rfr, (15) 
-oo 
f+oo 

Sfu{oo)= / e-''^^{f{t)u{t-r))dT. 



In practice, these noise terms are found not to vary much with frequency in the antenna's 
sensitive range. It is often convenient to parametrize these spectral densities with three 
values; noise temperature T^, noise resistance r„, and noise reactance x„. They are defined 



as 



T„ = -^5A-[Im(^;0P, (16) 
Kb 



Sf [lm{Sfu)Y 



(17) 



= (18) 

where ks is Boltzman's constant. For simplicity, we set lm{Sfu) = 0. Although, in general, 
the correlation between the force and velocity noise is non-zero, the effects of a non-zero 
Sfu can normally be accounted for by a renormalization of the transducer spring constant 
20| . The real part of Sfu is normally zero when a SQUID amplifier is used. 

Once a complete expression for a^uj) has been obtained, the energy SNR can be calculated 
from 

SNR = — / (r{uj)duj. (19) 

Note that we have consistently used a double-sided spectral density in contrast to the single- 
sided convention adopted by LIGO. By putting in numerical values for all parameters, this 
integral can be evaluated. For many of the parameters below, we chose values beyond what 
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has been demonstrated experimentally so as to represent an advanced spherical detector. 
Such an advanced detector could operate concurrently with the first LIGO interferometer as 
a result of aggressive research and development efforts now being planned ||T6[. Some of our 
parameter values are only slight extrapolations beyond currently demonstrated values, while 
others are instead upper bounds to the technologies being pursued. A detailed consideration 
of the research burden to meet each of our assumed values is beyond the scope of this paper. 
As stated in the introduction, our motivation is to clarify how such spherical detectors 
could complement the interferometer experiments by examining the detectability of high 
frequency events. We feel it is likely that resonant-mass detectors with an energy sensitivity 
approaching that derived with our parameters can be developed and built to operate within 
the time frame between the completion of the first LIGO interferometers and the operation 
of advanced interferometers. 

We chose to model aluminum spheres at a physical temperature of 50 mK, instrumented 
with six identical sets of three-mode inductive transducer systems located with the dodec- 
ahedral TIGA geometry WM. The lowest temperature that an aluminum bar antenna has 



been successfully cooled to is 95 mK Two-mode transducers are in use on a number 
of operating cylindrical resonant-mass antennas H-ffl and a three-mode system has been 



demonstrated at 4 K in a smaller, test antenna p6|. A constant mass ratio between the 
effective mass of the sphere and the intermediate mass as well as between the intermediate 
and transducer masses of 100:1 was used, and all mechanical quality factors (Qs) were as- 
sumed to be 40 X 10^. The highest mechanical Q that has been obtained in an inductive 
transducer is 24 x 10^ [ p^ . 



The transducer electronics were assumed to be a 9 cm diameter inductive pickup coil at- 
tached to a SQUID amplifier with a quantum-limited noise temperature, i.e., Tg = 1 huJo/ks- 



Quantum-limited SQUIDs have been constructed [^], but they are not useful for inductive 
transducers because of their low input coil inductance. Getting a suitable quantum-limited 
SQUID is an area of intense research. Wellstood's group at the University of Maryland is 
developing a quantum-limited SQUID for use in a gravitational wave transducer. The best 



noise temperature they have achieved in a SQUID with high enough inductance to couple to 
the transducer coils is Tg ~ 28 huj^/kB |2^. The prospect of approaching the quantum-limit 
in a practical SQUID in the next several years looks real. With proper matching, the trans- 



ducer noise is limited by the noise of the SQUID, so the value of T„ in Eq. (|T^) becomes 
equal to T^. The noise resistance is r„ = /c£;/47r/o, where fc^; is the real part of the spring 
constant fctrans that is due to the electrical interaction between the transducer mass and the 
pickup coils. The ratio kE/ktrans is the coupling between the electrical and mechanical parts 
of the transducer. For the value ks, we took the product 3.78 x 10® N/m^ x coil area, 
based on measurements made in our laboratory at Maryland . 



Taken together, these parameters define the overall sensitivity of the antenna. Energy 
sensitivity Eg is defined as 

Eg = (20) 

" SNR ^ ^ 

It is useful to express this sensitivity in relation to the standard quantum limit, the minimum 

sensitivity possible using a linear amplifier |31|]. Expressed as a multiple of this standard 



quantum limit, the antenna sensitivity becomes 



N=^. (21) 
nujQ 



As a comparison for the numerical result, we calculate an approximate value of from 

(22) 



N - — 

nujn 



T I 1 

+ - — — + — — + 



5 a2Qint asQtrans, 

Equation ( p2D is derived in the Appendix. The parameters 5 and are also defined and 
computed in the Appendix. Substituting the values of T, T^, Qcfr, Qint, and Qtrans assumed 
above into Eq. ( P^ gives 

A^ ^ 1.0 + 0.96 + 0.87 + 0.78 

^ 3.6. (23) 

We calculated SNRs for eight different spheres. The diameter of the lowest-frequency 
sphere was chosen to be the largest size that might be constructed, 3.25 m. The size of 
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the highest-frequency sphere was chosen so that its lowest quadrupole frequency coincides 
with the peak in the spectrum of the coalescing binary neutron star data. This peak is at 
twice the rotation frequency of the transient, barlike structure that forms immediately after 



coalescence This assumption gives a sphere diameter of about 1.05 m. The remaining 
sphere sizes were chosen to give reasonably continuous coverage of the frequency band 750 Hz 
to 2700 Hz. In addition to transducers tuned to the lowest quadrupole mode of the sphere, 
a system tuned to the first excited mode was examined. Coccia, Lobo, and Ortega [|l5l have 
shown that the cross section of the excited mode of a large sphere is 2.72 times greater than 
that for the lowest mode of a small sphere at the same resonance frequency. This allows the 
calculations of SNRs for the higher mode. 

Figures |l] and |^ show the sensitivities of the eight spheres in the ground state and 
excited state, respectively. These figures also show the sensitivity of the first LIGO and 
advanced LIGO interferometers for comparison. The data graphed is the strain spectrum of 
the detectors, h{u), defined as 



The strain spectrum is a measure of what frequency distribution an incoming gravitational 
wave would have to have in order to produce an output in a noiseless detector that mimics 
the output of the real detector's noise. It is a useful way to compare detectors because it 
is independent of source waveform and thus is solely a characteristic of the antenna. The 
energy SNR per unit bandwidth, a{uj), is related to the quantities he and /irms used by LIGO 
and defined in Ref. by 



Using these quantities, a value for the SNR of LIGO is often estimated as the maximum 
value of {he/ iV5hj.^s)Y , which gives a rough approximation to the integral in Eq. (^). The 
factor of is necessary to convert to "random direction and polarization" p3 |. 



These figures show that the spherical resonant-mass antennas have a sensitivity interme- 
diate between the first and advanced LIGO interferometers within a fractional bandwidth 
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of about 10% each. The collection of all the TIGAs, or the "xylophone", is a more sensitive 
detector than the first LIGO interferometers from 750 Hz to 2700 Hz in the lowest mode 
and from 1350 Hz to 5100 Hz in the first excited mode. In these frequency regions LIGO's 
sensitivity is constrained by photon shot noise in the lasers 0. 

III. SIGNALS AND RESULTS 

In order to integrate cr{uj) and obtain the SNR for the spherical antenna, numerical 
values for a gravitational waveform from an astronomical event are needed. For inspiraling 
and coalescing binary neutron stars, we used the waveform published by Zhuge, Centrella, 
and McMillan |^2|]. We Fourier transformed the time-domain data using the convention 
h{u!) = e~^'^^h(t)dt. Although the calculation by Zhuge et al. is among the most 
complete available, it is still only a first survey of how binary neutron stars may behave. In 
particular, it models gravity with a purely Newtonian formula. The inclusion of general- 
relativistic corrections may significantly change the orbit and lead to differences in the 
waveform. Our results derived using this waveform should be seen in this light. 

The frequency-domain waveform for the inspiral and coalescence phase of the binary 
neutron star evolution is shown in Fig. |[ Zhuge et al. generated the waveform using 
a three-dimensional numerical simulation which models the neutron stars as nonrotating 
polytropes. The neutron stars were chosen to have equal masses of 1.4 M© each, since all 
known cases of neutron stars in binary systems have this mass Q. Initially, the distance 
between the stars was chosen to be much larger than the diameter of individual stars, so 
tidal gravitational effects are negligible. Thus, the stars are originally spherical, with a 
radius of 10 km. The initial orbit was chosen to be nearly circular and it evolves due to 
Newtonian gravity with a frictional term added to simulate the energy loss to gravitational 
wave emission. When the stars spiral together, tidal distortions in each star's shape grow 
larger and the evolution approaches coalescence. 

Once the separation between the stars is comparable to the neutron star's radius, hydro- 



12 



dynamic effects become important and an approximation of the nuclear equation of state is 
required. Zhuge et al. used 

P = Kp^^^''' (26) 

as the equation of state, where P is pressure, p is density, is a constant that measures the 
specific entropy of the nuclear matter and n is the polytropic index. A value of n = 1 was 
used for the waveform we analyzed. Smooth particle hydrodynamics (SPH) is then used to 
model the coalescence phase once the equation of state is specified. 

The gravitational waveform was calculated from the complete orbit of the binary neutron 
star system using the quadrupole approximation by Zhuge et al. This approximation ignores 
contributions from mass moments higher than the quadrupole, but is valid for nearly Newto- 
nian sources . In the transverse-traceless (TT) gauge, the gravitational wave amplitude 
is 



"5" = ^'r. (27) 



where / is the second time derivative of the reduced quadrupole mass moment of the source. 
The amplitude of the "plus" and "cross" polarizations of the gravitational wave, expressed 
in spherical coordinates, are 

h+ = ~^{iee-iH) (28) 
= (29) 

The absolute scale of these amplitudes requires a choice for r, the distance from the detector 
to the source. We used r = 15 Mpc, the approximate distance to the Virgo cluster of galaxies 
33| . This distance is estimated by Phinney |p to be below the most optimistic value to get 



3 events per year, 23 Mpc. Scaling from this value, about 1 event per year is predicted at 
15 Mpc in this optimistic limit. 

The waveform from the Newtonian inspiral with friction was then meshed onto the 
waveform from SPH by Zhuge et al. to get a complete waveform for the whole binary 
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neutron star evolution. Since the orientation angles of the binary system are not known 
a priori, and in fact are values that the spherical antenna can determine experimentally 
p!8| , we averaged the waveform over these unknown angles. This averaging is done so that 
the energy per frequency, dE/df, radiated by the binary system is held constant. Thus, in 
Eq. (§), 

\h{co)\' = {\h+ico)\' + \h.ico)\'), (30) 

where (■ ■ ■) denotes an average over all source angles. It is this waveform that is shown in 
Fig. ^ and was used as input in Eq. (|14D. 

Once a numerical expression for the waveform h{uj) was made available to us, it was 
possible to obtain SNRs for the eight spheres with the diameters shown in Table |I[ To do 
this, the integral in Eq. ([191) must be evaluated. Performing this integration with the h{uj) 
from Eq. (^0]) gave the SNRs in the column marked "Total" in Table |. The row marked 
"Xylophone" is what an array of all eight TIGAs acting together could accomplish and is 
the sum of each SNR in the rows above. The row marked "First LIGO Interferometers" is 
for comparison with the first LIGO interferometers and was calculated by using the same 
waveform integrated with the strain spectrum published for LIGO Since the waveform 
is of finite extent in time, the frequency domain data is not accurate below 300 Hz. In 
order to get a reasonable value for the SNR of LIGO, it was necessary to extrapolate the 
data below this cut off and into LIGO's sensitive region. We did this with the analytical 



waveform in Eq. (44) of Ref. Note that this equation as published has a factor of 2 

error. The value n/12 should be instead 7r/6 [0]. We used the corrected version for the 
extrapolation. It was also necessary to divide the LIGO SNR by a factor of 5 to represent 
a wave with random direction and polarization |Q. This correction is unnecessary for the 
spherical antennas because they are equally sensitive in all directions. 

To determine how effective the TIGAs can be in observing the coalescence phase of the 
binary neutron star evolution, we separated the waveform into two parts. The inspiral occurs 
from t = s to t = 0.234 s, and the coalescence occurs from t = 0.234 s to t = 0.241 s. This 
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division of time was chosen so that the instantaneous frequency at t = 0.234 s coincides with 
/dyn, the dynamical instabihty frequency identified by Zhuge et al. as the frequency where 
the neutron stars cease to act as point masses. The separate time-domain waveforms were 
then multiphed by a Hahn windowing function I^S] before Fourier transforming, to ensure 
that the division was smooth and no spurious high-frequency signals were artificially created. 
The SNRs obtained from each of these separate waveforms are shown in the columns marked 
"Inspiral" and "Coalescence" , respectively, in Table |l[ 

With the results, we can calculate the energy sensitivity Es achieved by the spheres from 
Eq. (0). The energy deposited in the 3.25 m sphere in the lowest mode is 1.79 x 10^^^ J, 
from Eq. (1). Thus, with a SNR of 11.3, 

1.79 X 10-29J 



E, 



11.3 

3.0^a;o, (31) 



in good agreement with the approximate calculation in Eq. (p2D. 

We also calculated SNRs for the spheres instrumented with resonant-mass transducers 
tuned to the first excited mode of the antenna. The same physical parameters were used to 
model the spheres and the same waveforms used as signals. This data is shown in Table 0. 

Binary neutron star events are the best understood signals for gravitational wave detec- 
tion, but because of the high operating frequency of spherical antennas, other astronomical 
sources may be important for the spheres. One possible high-frequency signal is from the 
dynamical bar-mode instability of a rapidly rotating star. This event may be detectable 
by spherical antennas provided the star is compact enough. The stellar radius that this 
event occurs at is uncertain. This instability has been investigated by Smith, Houser, and 



Centrella and it is their numerical waveform data that we used. 

Figure ^ shows the spectrum for the bar-mode instability which develops following "cen- 
trifugal hangup" during the core collapse of a massive star. This gravitational waveform was 
generated by Smith et al. using three-dimensional numerical simulations which modeled the 
star polytrope, with equation of state 
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P = ^. (32) 

n 

Here P is pressure, p is density, e is specific internal energy and n is the polytropic index. A 
value of n = 3/2 was used for the waveform in our simulation. A total mass of 1.4 was 
assumed, as the star is expected to end up as a neutron star. An equatorial radius for the 
centrifugal hangup of 20 km was also assumed. This is the radius where the star's collapse 
is arrested by the centrifugal effects from the rotation. A realistic value for this radius is 
not known, but could range from as high as 3000 km down to as low as 10 km [^]. The 
initial rotation had a ratio of rotational kinetic energy to gravitational potential energy of 
r ~ 0.30. Newtonian gravity was assumed and the gravitational radiation produced from the 
dynamical bar-mode instability was calculated in the quadrupole approximation, as with 
the coalescence waveform. Back reaction from gravitational wave emission was ignored. The 



bar mode, i.e., m = 2 mode, was used as it is expected to be the fastest growing mode |^ 
This waveform was generated for the dynamical ha,T instability, which is driven by Newtonian 
hydrodynamics and gravity rather than the secular instability, which is due to dissipative 
processes such as gravitational radiation reaction. The dynamical instability develops on 
a time scale of about one rotation period while the secular instability grows over several 
periods or even more slowly. The choice of r = 0.30 is just above the dynamical stability 
limit of 0.2738 and thus is a reasonable approximation for a star that spins up, due 
to collapse or accretion, and becomes dynamically unstable The star's evolution was 
simulated by Smith et al. using SPH and from this evolution the gravitational waveform is 
calculated by using Eq. (|27|) . 

The choice of r, the distance from source to detector, is not as simple as for the binary 
neutron star. There is less observational evidence for stars with bar-mode instabilities. Such 
rapidly rotating stars may be formed from supernovae, so the rate of supernovae might be 
taken as a reasonable guide to the rate of this gravitational wave event. Out to the Virgo 



cluster of galaxies (15 Mpc), the supernovae rate is estimated at a few per year ||3^. We 



took the source distance to be r = 1 Mpc, which has an estimated event rate of 1 every 10 
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years |3^. Once the magnitude of each polarization state was evaluated, the same average 



over angles as in Eq. ( pOD was performed to give the waveform shown in Fig. ^. The SNRs 
were calculated by using the method described in Section 0. These values are shown in 
Table ^ for both the ground state and the first excited quadrupole modes. The 1.45 m, 
1.25 m, and 1.05 m diameter spheres do not have data listed for the excited mode because the 
frequency-domain waveform cuts off at 3500 Hz, which is below the resonance frequencies 
of these spheres. This is because the granularity of the time domain data provided was 
too great for frequencies above 3500 Hz. However, we believe it is safe to assume that the 
frequency-domain data would remain below h{uj) = 10~^^ Hz^^ and thus the SNRs for these 
two spheres would be negligible. 



The secular instability can develop for r values greater than 0.1375 After the 

dynamical instability has run its course, Houser et al. find the system has evolved to a 
nearly axisymmetric state with a core having r = 0.26, which is above the secular instability 



limit. Lai and Shapiro |^ have investigated the secular instability which may develop in 
which a Maclaurin spheroid evolves into a Dedekind ellipsoid producing a gravitational wave 
signal sweeping in frequency from possibly near 1 kHz down to zero. A different type of 
secular evolution would apply if the calculations of Durisen et al. or Williams and Tohline 



4l| correctly predict the end point of the dynamical instability as a bar surrounded by a 



ring, rather than the spheroid found by Houser et al. If a rapidly spinning bar is produced, 
the secular evolution changes the bar from a Jacobi-like ellipsoid into a Maclaurin spheroid 
with a gravitational wave signal from 500 Hz to as high as 3 kHz. 

Lai and Shapiro give analytical waveforms for these two different secular instabilities. For 
the same I.AMq star at 1 Mpc with a radius of 10 km averaged over all source and detector 
angles, the Dedekind waveform would give an energy SNR of 2000 in the 3.25 m diameter 
sphere and 1200 for the first LIGO interferometers. Such a strong signal would be detectable 
even at the Virgo cluster distance. However, it seems likely that the starting frequency for 
the Dedekind evolution will be well below the spherical detector's band unless the star is 
initially spinning just below the dynamical instability limit and is nearly incompressible. 
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In order for the starting frequency of this event to exceed the 795 Hz frequency of the 
3.25 m sphere, Lai and Shapiro's Fig. (5) indicates one needs a polytropic index n < 0.7 for 
r = 0.26, while even n = does not suffice for r < 0.24. In contrast, the frequency span of 
the signal from Jacobi-like evolution seems to have a higher possibility of overlap with the 
band covered by the spheres; however, the chief uncertainty for that signal is whether or 
not the dynamical instability produces a spinning bar rather than a spheroid. // Jacobi-like 
evolution occurs and its starting frequency is below 800 Hz, then the signal computed by 
Lai and Shapiro is a very strong source for a big sphere: we calculate an energy SNR of 4000 
for a 3.25 m diameter sphere and an energy SNR of 240 for the first LIGO interferometers 
from a waveform averaged over all angles and a source distance of 1 Mpc. 



IV. CONCLUSION 

The results in Table | for the spherical antenna tuned to the lowest quadrupole frequency 
interacting with gravitational radiation from binary neutron stars shows that spherical an- 
tennas operate at a level that is complementary with the first LIGO interferometers. The 
largest sphere obtains an energy SNR of 11.3 at a distance of 15 Mpc. If a SNR of 10 in 
this detector is a sufficient threshold for a three-way coincidence experiment, then a source 
with angle averaged strength could be detected out to a distance of 15.9 Mpc. With the 
most optimistic estimate of the coalescence rate, 3 events per year out to 23 Mpc 0, the 
expected rate of detection is 

1 /l5.9Mpc\^ , , 

Td = 1.15 X 3 jT-^ (33) 

^ \ 23 Mpc J ^ ' 

= 1.15yr-^ (34) 



The factor of 1.15 is due to statistical preference for angles that give high SNRs (see ||33 
for details). For detection of gravitational radiation from binary neutron stars at a distance 
of 15 Mpc, a 3.25 m diameter aluminum sphere near the standard quantum limit will be 
sufficient. The upper limit on the event rate at this distance is about 1 coalescence per year 



Table | also shows that a large sphere instrumented at the lowest quadrupole frequency 
does not hold out much hope of seeing the details of binary neutron star coalescence. Even 
the 1.05 m diameter sphere, whose size was chosen so that the lowest quadrupole mode was 
at the maximum of the coalescence wave spectrum, does not manage to reach a SNR of 1. 
As the frequency of the sphere goes up, the radius, and with it its mass, goes down. At 
frequencies where the waveform from Zhuge et al. is strong, the energy cross section of the 
sphere is too small to detect much. This raises the question of the reliability of the numerical 
waveform, especially of the frequency /peak associated with the barlike transient. According 
to Centrella |^2|, the qualitative shape of the waveform is fairly reliable, but the exact 



position of this peak and other structures may change as numerical relativity techniques 
improve. If /peak were to be found at a lower frequency, closer to the lowest quadrupole 
mode of one of the larger spheres, the prospect for a SNR greater than 1 for the coalescence 
phase might improve. 

The data in Table ^ for spheres sensitive at the first excited quadrupole mode to inspi- 
raling and coalescing binary neutron stars appears a little more promising for the detection 
of coalescence. The largest sphere still has the highest overall SNR, but it is reduced from 
its value in the ground state. The 2.00 m diameter sphere, with an excited quadrupole 
frequency of 2483 Hz, is now the sphere tuned to /peak- It does not quite manage a SNR of 
1 either, but it does have a higher SNR for the coalescence than for the inspiral phase, as 
do the 1.70 m and 1.45 m spheres. 

The data in Tables | and taken together, suggest that a xylophone of spheres acting 
collectively and in collaboration with LIGO may be the best approach to detection of binary 
neutron star coalescence. The smaller spheres do not contribute much to the xylophone SNR, 
so four spheres, from 3.25 m to 2.00 m in diameter, would be enough to give most of the 
xylophone benefits. If these four spheres could be instrumented at both the ground and 
the first excited quadrupole frequencies, a fairly wide spectrum, from 750 Hz to 2700 Hz 
continuously, could be monitored. The spheres could rely on LIGO to provide a high SNR 
detection of the inspiral. Once the event is discovered, the spheres would have the less 
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strenuous task of identifying and analyzing the high frequency coalescence. Thus, even the 
modest SNRs found in our work may still prove useful in gathering astrophysical data on 
coalescing binary neutron stars. 

The data for the rapidly rotating star shown in Table ^ is more encouraging. For the 
ground state modes, the 1.45 m sphere has its frequency near the peak of the spectrum and 
obtains a strong SNR of almost 10. The two spheres with quadrupole frequencies below the 
resonance of the 1.45 m sphere also have strong SNRs, about 7 and 5. Further away from 
the peak, SNRs fall off rapidly, especially on the high-frequency end. The SNR goes from 
0.935 for a sphere diameter of 1.25 m to well below 1 for a diameter of 1.05 m. Thus with 
only two antennas the peak signal can be easily found, provided their sensitive frequencies 
occur at the appropriate positions. 

The first excited state data is similar to the ground state, showing strong SNRs when the 
sphere's quadrupole frequency is at or near the 1765 Hz peak. In the excited mode, however, 
this occurs between the 3.25 m and the 2.75 m diameter spheres. With larger masses, these 
spheres have higher energy cross sections and thus obtain much higher SNRs. A SNR of 
25, from the 2.75 m sphere in the excited state, represents such a strong signal that the 
source position on the sky could be located to within almost 0.13 steradian |T^. As with 



the ground state data, there is a sharp drop in SNR to effectively zero about 500 Hz above 
the peak frequency, making location of the peak frequency possible. 

The rapidly rotating star waveform was generated by assuming a total mass of 1.4 Mq 
and a centrifugal hangup at 20 km. The location of the peak frequency, which is twice the 
rotation frequency of the star, is very sensitive to the values of these parameters. It can be 



as low as 1 Hz for a 1.0 M0, 3000 km star up to 6000 Hz for a 2.0 M©, 10 km star 
Since the appropriate values for these star parameters are not known, and in fact are values 
we could hope to determine from gravitational wave data, actual signals from this source 
could potentially be outside the sensitive range of spherical antennas. A peak frequency 
above about 2500 Hz, corresponding to a 1.4 Mq, 10 km star, would be extremely difficult 
to detect outside our galaxy. This would limit the number of events to a few per century. 
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Detecting these higher frequency signals depends on the accuracy of the current template, 
especially the secondary peak in Fig. ^at 400 Hz. Many details of the rapidly rotating star's 
evolution are not well understood and this waveform may undergo substantial changes as 
the field of numerical relativity advances. 

A xylophone of a 3.25 m, a 2.75 m, a 2.35 m and a 2.00 m diameter sphere instrumented 
at both the ground state and the first excited quadrupole state working along with LIGO, as 
suggested for the inspiraling and coalescing binary neutron stars, would do an effective job 
of searching for rapidly rotating star events. For favorable mass and radius parameters, the 
excited state of the 2.75 m diameter sphere would be sensitive, with a SNR threshold of 10, 
out to 1.6 Mpc. This will be sufficient to observe one event every decade, provided nearly 



all type lb and type II supernovae are sources of this event |^. This xylophone would also 
determine the star's rotation frequency, as a large SNR would be seen in the more massive 
sphere and effectively nothing would be seen in the smaller. This would locate the peak 
frequency, and hence the rotation frequency, to within a few hundred hertz. Negative results 
from such a xylophone would restrict the parameter space, providing data about neutron star 
development and equation of state. The xylophone covers a frequency range corresponding to 
hangup radii between 45 km down to 15 km. The broadband interferometers have difficulty 
detecting this dynamical instability, but do have a high SNR for the secular instability. Since 
these two gravitational wave events may occur in succession, a collaborative effort between 
spheres and interferometers would prove effective as with the binary neutron stars. 

In this paper, we have compared broadband interferometers with resonant-mass anten- 
nas for detection of high frequency gravitational radiation. Another possible technique for 
detecting high frequency events involves dual-recycled interferometers ||4^. This approach 



allows laser interferometers to become much more sensitive within a narrow bandwidth at 



the expense of sensitivity elsewhere. Krolak, Lobo, and Meers have looked at SNRs for 
inspiraling neutron stars interacting with recycled interferometers using a simplified strain 
spectrum and an analytical formula for the waveform. Their results for the unrecycled 
LIGO are in close agreement with ours. Further investigations of SNRs comparing spherical 
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antennas with dual-recycled interferometers is being done at Caltecli p4 . 

In addition to the astronomical sources of gravitational radiation that we investigated 
here, there may be other high frequency signals potentially detectable in detectors optimized 
for the 1-2 kHz band. Events that might produce high-frequency gravitational waves in- 
clude coalescence of a neutron star with a black hole or a black hole with a second black 
hole asymmetric core collapse and bounce in supernovae spinning neutron stars 



4J], and cosmic string vibrations Especially promising may be the black hole coales- 



cences and spinning neutron stars. Excitation of the high-frequency (f ~ 1200 — 3250 Hz 
(IOM0/M) black hole quasi-normal modes would give a relatively strong signal at kilo- 
hertz frequencies. A gravitational-wave antenna detection of this radiation could provide 
the observational "smoking gun" to confirm the existence of black holes. Experimental evi- 
dence of gravitational radiation from black hole coalescence would undoubtedly also provide 
great insight into relativistic gravity. Spinning neutron stars are a periodic source that could 



radiate strongly for months The frequency of the waves would be twice the rotation 
period, often above a kilohertz, and details of the wave would tell much about the structure 
of neutron stars. According to Thorne ^M, "the deepest searches for these nearly periodic 
waves will be performed by narrow-band detectors ... e.g., dual recycled interferometers 
or resonant-mass antennas." We call on the numerical relativity community to continue 
to develop reliable waveforms for all possible high-frequency events. It is only through the 
combined efforts of everyone; interferometer and resonant-mass experimentalists, as well as 
numerical and analytical theorists, that confirmed, direct detection of gravitational radiation 
will become a reality. 
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APPENDIX: ASYMPTOTIC EXPANSION OF THE ENERGY SENSITIVITY 



Equation (pO[) defines the energy sensitivity Eg as the ratio of E, the energy the signal 
would deposit in a bare antenna initially at rest, to the SNR calculated by integrating (j{uj) 



given by Eq. (14). Although we calculated Eg directly in this fashion, useful insight into 
the dependence of Eg on the detector parameters comes from considering the asymptotic 
behavior for the Qi approaching infinity as indicated by Eq. 

First, consider the lossless case when all Qs are infinite. Then Re(?/22) is zero, and Price 
2I| shows that the expression for (j{uj) can be rewritten exactly as 



a{uj) = (Al) 



where 



enico) = I 'f^'^'' r„, (A2) 
Tn is the mechanical amplifier noise temperature defined by Eq. (|16D, and 

Z'li — '^n ~\~ jXn (-^3) 

is the mechanical amplifier's complex noise impedance, defined by Eqs. (p!7|JT8| ). 

As Price explains, e„(to') has the following physical interpretation: it is the spectrum 



of the energy which would be dissipated in r„ if the signal were applied when the amplifier 
had been replaced by its input impedance plus an additional mechanical impedance equal to 
Zn- If fi{uj) varies little over the band where e„(ti;) is large, then the signal is approximately 
an impulsive force which deposits the energy 

E = (A4) 

initially all in antenna motion. Since the hypothetical detector with z„ inserted is assumed 
lossless except for r„, eventually 100% of the energy E ends up being dissipated in r„. Thus 
from this physical argument we have 
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— e„(^) ^ — , (A5 

independent of the values of r„ and the detector masses and springs. Therefore, in the 
lossless case, the SNR is 



coo (^QJ 

SNR= / —(t{uj) (A6) 
-oo zvr 

- dij enjiv) ^^^^ 

-oo 271 ksTn 
KB^n 

provided /i and T„ vary little over the frequencies where e„ is significant. Hence, for an 
impulsive force. Eg = fc^Tn in the lossless case independent of the other detector parameters. 

Next, consider the case where all Qs are infinite except that for the antenna. Then 
the antenna dissipation produces a thermal force noise acting on rries with spectral density 
2kBTmcsuJo/Qc{i- Therefore, Eq. ([T^) becomes 

''^'^> - Su + Sf\y22\' + 2Re{y22Sfu) + 2kBT^\y,,r ^ ' 

For convenience, we define the energy sensitivity expressed in temperature units to be Tp, 
the detector's pulse-detection noise temperature: 

keTp = Es. (AlO) 



Motivated by Eq. (26) in Price we use Eqs. ([A1| , [A2|) to rewrite Eq. ([A9|) in terms of 



e„(ti;), and then expand the integral for SNR to first order in Q^^ and use Eq. (|A5|): 

SNR = (All) 



. 2n kBT„, 1 + J^o ^"(") 
E ( Tluq \ 



ksTn \ TnQes25fo 



(A13) 



where 



2<^/0 = .oo , xn2 (^14) 
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defines the fractional bandwidtli S of the lossless detector in a way which arises naturally in 
the SNR expansion, and which agrees with the intuitive definition of bandwidth if en{uj) is 
box-shaped. 

Thus when all other Qs are infinite, the asymptotic form of Eg for large Qes is described 

by 

1^1/ TUq \ ^^^^^ 

Tp Tn \ TnQcs'25fo J 
T 

Tp^T„ + --^7r asgcff^oo. (A16) 

When other Qs besides Qes are also finite, analysis along the lines above is more cumbersome, 
but one finds 

rp N 1 

d ^ aiQi 



where the constants ai can be determined numerically. Comparison with Eq. ( |A16|) shows 



that ai = I/tt always. Thus, a higher fractional bandwidth 6 reduces the effect of finite Qs 
T 

Price |PD| argues that a clever choice of the masses and springs to maximize 5 for a 



on Tp. 



given value of r„ is the optimal design strategy, and that a design in which e„(co') is made 
"optimally flat" is nearly optimal in this sense. However, in this paper, we have chosen the 



simpler design originally proposed by Richard which has constant successive mass ratios 



since we find both designs give similar values for Tp. 

The values of a, for i > 2 depend on the masses and springs chosen for the detector 



design. The values for ai,a2, and as that we calculated and used are shown in Table |V. 
These values are for the case of a constant mass ratio between stages of 100. This ratio 
was chosen as a convenient value that gives a final transducer mass slightly higher than the 
mass that gives the lowest Tp. A higher mass is more resistant to uncertainties in noise 
resistance and thus is preferred experimentally to a lower than optimal mass. The design 
rule suggested by Richard for choosing a constant mass ratio given an r„ and nies can 
be written as 
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2iV-l 

^] ' =C^^, (A18) 

where iV is the number of modes of the system and C is 2. The design we chose corresponds 
to = 3 and C = 2.61; the in Table |V| apply for any such design. 
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TABLES 

TABLE I. Energy signal-to-noise ratios for binary neutron star evolution in the lowest ^ = 2 
mode of the sphere. The distance to the neutron stars is taken to be 15 Mpc and each neutron star 
has a mass of 1.4 M0 and a radius of 10 km. The waveforms have been averaged over source and 
detector angles. The individual spheres have a sensitivity about 3 times the standard quantum 
limit. The row labeled "Xylophone" is obtained by summing the signal-to-noise ratios for each 
sphere. 



Diameter 


Frequency 


Coalescence 


Inspiral 


Total 


3.25 


795 Hz 


0.0113 


10.6 


11.3 


2.75 


940 Hz 


0.00985 


5.79 


6.28 


2.35 


1100 Hz 


0.0146 


3.43 


3.88 


2.00 


1292 Hz 


0.00948 


1.40 


1.64 


1.70 


1520 Hz 


0.00853 


0.907 


1.09 


1.45 


1782 Hz 


0.0104 


0.558 


0.719 


1.25 


2096 Hz 


0.0197 


0.285 


0.449 


1.05 


2461 Hz 


0.126 


0.0886 


0.407 


Xylophone 




0.210 


22.8 


25.6 


First LIGO Interferometers 


0.00406 


58.2 


58.8 
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TABLE II. Energy signal-to-noise ratios for binary neutron star evolution in the first excited 
1 = 2 mode of the sphere. The waveforms have been averaged over source and detector angles. The 
distance to the neutron stars is taken to be 15 Mpc and each neutron star has a mass of 1.4 and 
a radius of 10 km. The individual spheres have a sensitivity about 3 times the standard quantum 
limit. The row labeled "Xylophone" is obtained by summing the signal-to-noise ratios for each 
sphere. 



Diameter 


Frequency 


Coalescence 


Inspiral 


Total 


3.25 


1528 Hz 


0.0235 


2.43 


2.93 


2.75 


1806 Hz 


0.0289 


1.45 


1.88 


2.35 


2113 Hz 


0.0593 


0.679 


1.13 


2.00 


2483 Hz 


0.349 


0.222 


1.09 


1.70 


2921 Hz 


0.0688 


0.0448 


0.224 


1.45 


3425 Hz 


0.00709 


0.00591 


0.0252 


1.25 


3973 Hz 


0.00111 


0.00150 


0.00513 


1.05 


4729 Hz 


0.000674 


0.000733 


0.00279 


Xylophone 

First LIGO Interferometers 


0.538 
0.00406 


4.83 
58.2 


7.29 
58.8 
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TABLE III. Energy signal-to-noise ratios for the rapidly rotating star waveform for ^ = 2 modes 
of the sphere. The waveforms have been averaged over source and detector angles. The dashes for 
the excited mode of the 1.45 m, 1.25 m and 1.05 m spheres represent the fact that the signal data 
cuts off below the resonance frequency of these modes. The distance to the star is taken to be 1 
Mpc, the mass of the star is taken to be 1.4 Mq, and the radius of centrifugal hangup is taken to 
be 20 km. The individual spheres have a sensitivity about 3 times the standard quantum limit. 
The row labeled "Xylophone" is obtained by summing the signal-to-noise ratios for each sphere. 
Diameter Ground State Excited State 





Frequency 


SNR 


Frequency 


SNR 


3.25 


795 Hz 


0.0661 


1528 Hz 


19.5 


2.75 


940 Hz 


0.220 


1806 Hz 


25.0 


2.35 


1100 Hz 


1.22 


2113 Hz 


1.08 


2.00 


1292 Hz 


4.75 


2483 Hz 


0.00434 


1.70 


1520 Hz 


6.95 


2921 Hz 


0.00278 


1.45 


1782 Hz 


9.91 


3425 Hz 




1.25 


2096 Hz 


0.935 


3973 Hz 




1.05 


2461 Hz 


0.00168 


4729 Hz 




Xylophone 




24.1 




45.6 



First LIGO Intorfcromotors 0.197 0.197 
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TABLE IV. Coefficients in the asymptotic expansion of energy sensitivity. The Oj determine 
the influence of the Qs on the noise temperature. The value S is the fractional bandwidth of the 

detector in the lossless limit. 



Constant Calculated Theoretical 
ai 0.34 I/tt = 0.32 

a2 0.41 
a3 0.40 

S 0.103 A/mtrans /^"■int = . 1 00 



31 



FIGURES 
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FIG. 1. The strain spectrum of the eight spherical antennas in the lowest quadrupole mode, 
shown with solid and dashed lines. The different line styles have no significance other than to 
differentiate the separate strain spectra. The upper dotted line shows the strain spectrum for the 
first LIGO interferometer and the lower dotted line shows the strain spectrum for the advanced 
LIGO interferometer, for reference. The spherical antennas, each with a sensitivty about 3 times 
the standard quantum limit, are more sensitive than the first LIGO interferometers in a bandwidth 
of about 100 Hz to 300 Hz each and together span a total bandwidth from 750 Hz to 2700 Hz. In 
this band, the spherical antennas are a little less sensitive than the advanced LIGO interferometers. 
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FIG. 2. The strain spectrum of the eight spherical antennas in the first excited quadrupole 
mode, shown with sohd and dashed fines. The different fine styles have no significance other than 
to differentiate the separate strain spectra. The upper dotted line shows the strain spectrum for 
the first LIGO interferometer and the lower dotted line shows the strain spectrum for the advanced 
LIGO interferometer, for reference. The spherical antennas, each with a sensitivity about 3 times 
the standard quantum limit, are more sensitive than the first LIGO interferometers in a bandwidth 
of about 200 Hz to 600 Hz each and together span a total bandwidth from 1350 Hz to 5100 Hz. 
In this band, the spherical antennas are about equal to the sensitivity of the advanced LIGO 
interferometers. 
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f (Hz) 

FIG. 3. The frequency-domain gravitational waveform averaged over source orientation from 
the inspiral and coalescence of two neutron stars with mass 1.4 Mq and radius 10 km at 15 Mpc 
from the antenna. The sharp cut-off at 300 Hz is due to the finite extent of the time-domain data. 
The spectrum from 300 Hz to about 1000 Hz is due mainly to the inspiral phase. The frequency 
/dyn = 1566 is the dynamical instability frequency. The peak at / ~ 2500 Hz is associated with 
the transient barlike structure that forms immediately following the onset of coalescence |p2| . 
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10 




f (Hz) 

FIG. 4. The frequency-domain gravitational waveform averaged over source orientation from 
the bar-mode instabihty of a rapidly rotating star of mass 1.4 Mq and radius at centrifugal hangup 
of 20 km at a distance of 1 Mpc from the antenna. The primary peak at 1765 Hz is twice the 
rotational frequency of the star. 
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